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Abstract 

R-parity and L violation in the MSSM would be the origin of the neutrino 
oscillation observed in Super-Kamiokande. A distinctive feature of this frame- 
work is that it can be tested in colliders by observing decay products of the 
destabilized LSP. We examine all the possible decay processes of the neu- 
tralino LSP assuming the bilinear contribution to neutrino masses dominates 
over the trilinear one which gives rise to the solar neutrino mass. We find 
that it is possible to probe neutrino oscillations through colliders in most of 
the R-parity conserving MSSM parameter space. 
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I. INTRODUCTION 



Lepton (L) or baryon number (B) conservation is not a consequence of gauge invariance 
in supersymmetric extension of the standard model. The presence of both L and B violation 
would lead to fast proton decay, and thus one usually imposes a discrete symmetry which 
ensures L and/or B conservation. A typical example is R-parity discarding renormalizable 
L and B violating operators, as a consequence of which the lightest supersymmetric particle 
(LSP) is stable and a good candidate of dark matter. Another interesting possibility is 
to impose baryon parity (B-parity) which allows the following L violating terms in the 
superpotential of the minimal supersymmetric standard model (MSSM): 

l^iHuLi + —XijkLiLjE'^ + X'^jf^LiQjDl . (1) 

For the later use, it is convenient to generalize the bilinear term to include the usual /i-term: 
liaHuLa where a = 0, 3 and Lq = Ha- It remains an open question what kind of discrete 
symmetry can be derived from a more fundamental theory beyond the MSSM. In view of 
gauge symmetry origin of discrete symmetries, it has been shown that the standard R-parity 
and B-parity can only be consistent with the MSSM and the latter is preferred to arrange 
for proton stability without fine-tuning taking into account higher dimensional operators [|l| . 

R-parity (and L) violation in the MSSM would be the origin of neutrino masses and 
mixing, in particular, the recently observed atmospheric neutrino oscillation 0]. In the 
presence of the |AL| = 1 operators in the superpotential (|ID as well as in the soft terms, 
neutrinos get nonzero masses through the well-known mechanism of particle-sparticle mixing 
and one- loop radiative correction 0. A distinctive feature of this mechanism is that the 
LSP is destabilized and its decay modes can be easily identified in colliders. The most 
characteristic events for those L-violating phenomena are same-charge final-state dileptons 
due to the Majorana properties of two produced LSP's, which have little Standard Model 
backgrounds 0. Therefore, detecting the LSP decay, one can directly measure the quantities 
related to the L violating parameters f^J^, which might enable us to extract some crucial 
information on the neutrino sector of the theory and test neutrino oscillations in colliders. 

The purpose of this paper is to provide a detailed analysis of the decay processes of 
the LSP (preferably taken as the lightest neutralino) through both the bilinear and the 
trilinear L violating couplings in the MSSM whose sizes are determined by current neutrino 
experiments. The bilinear terms play an important role in the whole discussion. As they 
are typically the dominant sources for neutrino masses, our analysis will be based on the 
assumption that the atmospheric neutrino masses and mixing are mainly due to the bilinear 
terms. The collider signatures of the bilinear terms have been studied in Refs [[^-[T0[| taking 
only one lepton generation and the corresponding (tau) neutrino with a large mass ~ MeV. 
The correlation between collider signals and the atmospheric neutrino oscillation has been 
first analyzed in Ref . |^ and also in . In this paper, we extend these works to include fully 
the three lepton generations and analyze both two and three body decays of the LSP. The 
latter become important when the trilinear couplings are introduced to explain the solar 
neutrino oscillation. 

In the next section, we provide a comprehensive formulation of two and three body decay 
rates of neutralinos through both the bilinear and trilinear couplings. In section III, we 
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estimate the sizes of the R-parity and L violating parameters for which the current neutrino 
data can be accommodated. Then, we examine in section IV the region of neutrahno mass 
parameters for which the LSP decay and its branching modes can be measured. An emphasis 
is put on how to separate out the information on the bihnear terms from the contribution of 
the trihnear couphngs which are present to give the masses and mixing for the solar neutrino 
oscillation. Conclusion is given in section V. 



II. PARTICLE-SPARTICLE MIXING AND LSP DECAYS 

One usually rotates away the term HuLi in Eq. so that L violation only appears in the 
trilinear couplings whose consequences in collider experiments are studied extensively [11 . 
But, there also exist soft supersymmetry breaking bilinear terms which cannot be rotated 
away in general. These terms may play more important roles for the LSP decay P,[7|-p!0 . 
The relevant soft terms in the scalar potential are 

Vsoit = mjfjHdl^ + ml^m^ + {BHuHd + BiHuU + ml^H^L.Hl + h.c.} + ■ ■ • , (2) 

from which it is clear that the L violating soft parameters Bi,m\^^^ cannot be rotated 
away unless the soft parameters satisfy certain relations called the alignment condition Wl . 
This kind of alignment condition is usually achieved at the mediation scale of supersymme- 
try breaking which is presumed to be L conserving. However, at the weak scale where the 
minimization of the scalar potential takes place, the alignment breaks down by the renormal- 
ization group evolution of the soft parameters. Therefore, unless one enforces the alignment 
at the weak scale by some mechanism, the L violating terms being linear in Lj induce nonzero 
vacuum expectation values (VEVs) of sneutrino fields which are non-aligned with /ij. As 
a result of this misalignment, the mixing between particles and sparticles, in particular, 
neutrinos and neutralinos, arises to lead important phenomenological consequences. 

In the generic presence of the bilinear L violating terms /i, and nonzero sneutrino VEVs, 
the neutralino-neutrino mass matrix in the basis {B, W3, Va) is given by 



M 



N 



/Ml MzSwVu/v -Mzswvjv\ 

Ml -MzCwVu/v MzCwVa/v 

-yUa 

V / 



(3) 



Here, Mn is a 7x7 symmetric matrix, Vu = {Hu), vq = (Hd), Vi = (z>j) with i = 1..3, 
and sw, cw are the weak mixing elements. The neutralino-neutrino mass matrix can be 
diagonalized by the rotation A^*MnA^^. That is, the index i of the diagonalization matrix Nij 
runs for the mass eigenstates, and the index j runs for the weak eigenstates {B, W3, H^, z/^). 
We would like to remind that for the mass eigenstates Xi", the states with i = 1, 2, 3 are the 
neutrinos and i = 4 denotes the lightest neutralino. 

The char gino- charged lepton mass matrix in the basis (W~^ , H:^ , E^) and {W~,l~) is 

/ M2 V2MwVa/i 

Mc = V2MwVu/v fia I , (4) 

V M' 
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where is the 3x4 matrix with the components; M^q = —m{vk/vo, M^^ = 5ikm\. Two 
diagonahzation matrices Ur^l are defined by UrMcUI where Ur {Ul) casts the positively 
(negatively) charged states into the mass eigenstates. The first index i of the matrix Unj 
and U^ij run for the mass eigenstates; the states with i = 1,2,3 are the charged leptons 
(e, yU, r) and i = 4 and 5 the two charginos. 

In terms of the mass eigenstates, the neutral and charged current interactions take the 
form: 



(5) 



where 



^iI = Ee<xiV-^<x with e. = (l,-l,....,-l) 



cr=3 



(6) 



To denote the mass eigenstates in a conventional way, let us introduce the following nota- 
tions: Vi = x'i with i = 1 — 3 and Xi = Xi+s for i = 1 — 4. 

From Eq. (j^), one gets the following rate of the decay process x° - 
the neutralino x° is heavier than the Z and/or W boson: 



iZ or IfW^ when 



^/ 



X 

Mir 

Ml, 



qN 



+ 



(7) 



where I{x) = (1 — a;)^(l + 2x). Note that we have neglected the possible neutralino decay 
into the lightest Higgs whose rate is suppressed compared to the decay into W boson by 
factor of the order of M^/m|o- There could exist other LSP decay modes such as x° ^1 
which occurs at the one-loop level and has the branching ratio at most a few percent [|T^ . 

For the sake of discussion it is useful to use a seesaw formula valid for t>j/t>o, ^ 1. 
This allows us to decompose the diagonahzation matrices into two parts; the usual neutralino 
and chargino mixing, and the mixing between neutrinos (charged leptons) and neutralino 
(charginos) which arise from L violation. In this decomposition, the elements for the 
particle-sparticle mixing in Eq. (O) take the following factorized forms H : 



nN 



N 



.a=l 



(8) 
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where = Vi/vq — fii/fi. Note that no mixing effect can arise if the sneutrino VEVs are 
ahgned with /Xj, that is, = as discussed in the beginning of this section. The coefficients 
c^, and are given by [|6|,p!4| 

Ms Ml Mz Mz\ 




M^' M^' ^ /i ' 



V2,2s,^] , (9) 



/ ^,M2 1 , 2 , 

K2(— + -), — ( 



M| ^ ^^Mw. 



where Fn = MiMs/M;^ - Mfss/j/yU, Fc = Ms - M^S2/3//i, and M^ = c^Mi + s^Mg. In 
Eq. , and Cs have extra /ij / /i contributions fl^ which are irrelevant since they are 
canceled out in the physical process of neutralino decays as can be seen from the expressions 
in Eq. (|). Eqs. (|]) and (P) provide clear understanding of a few properties. The ratios 
^i/^j can be determined by measuring the branching fractions of the decay modes x° hW 
independently of the other (L conserving) supersymmetric parameters like M2 or /i also 
independently of the sizes of the trilinear couplings. The rate r(x° IW) depends on the 
charged lepton masses through (^. This mass dependence is suppressed by m\/Fc, but 
can be significant for Higgsino-like LSP and for large tan/3 as cf^ is proportional to l/c/3. 

When the LSP is lighter than W boson, it has only three body decay modes which have 
two major contributions: one is the bilinear contribution coming from the mixing Ur^l 
through the exchange of W or Z boson and the other comes from the trilinear L violating 
couplings A, A' through exchanges of squarks or sleptons. For the latter contributions, there 
exists a complete formulation in the literature [0. There are also contributions due to the 
charged and neutral Higgs exchanges which is neglected in our computation since they are 
generically smaller than the corresponding contributions from the W, Z exchanges. More 
important one comes from the slepton-Higgs mixing |]7|,p!6[| whose size is determined by the 



bilinear parameters fii, Bi and ^'l-h ■ This mixing effect can be translated to the trilinear 
couplings in the mass basis, which is included in our analysis. 

In order to calculate the rates of the neutralino decay x° — fifjfk, we recast the decay 
amplitudes in terms of the 12 matrix elements into the general form: 

Mix' ^ Mjfk) = Sjj{s,t,u)[v^Pm][ukPjVj] (10) 

+ Tij{s, t, u) [v^a^^PiVi] [uk(J^yPjVj\ , 

where I, J = L, R, and s = {p^ — piY, t = {p^ — pj)^ and u = {p^ — PkY- In our case, the 
neutralino decay has the modes; x° ~^ vvv., I'll, vqq and Iqq' . For each mode, nonvanishing 
matrix elements are presented in Appendix. The modes x^ —>■ vvv and vuu come solely from 
the bilinear contribution, and thus has one and two matrix elements, respectively. Note that 
the bilinear contributions appear in Slr or V/j, and only Sll and T^^, terms can interfere 
with each other neglecting the final state masses in our cases. It turns out this interference 
is negligible in the cases we will consider. 



5 



The mode x'^ — * ull contains 9 flavor modes of the charged leptons. Thus one could 
get information on the flavor structure of the trilinear couplings Xijk by identifying dilepton 
signals with missing energy. In a detector, the modes x° ~^ ^QQ and Iqq' can be seen as 
two jets with missing energy [uj]) and charged lepton with two jets (Ijj)- In the case of 
three body decays of the LSP, the values could be measured by separating out the bilinear 
contribution in the kjj decay modes. But the three body LSP decay length from the bilinear 
contribution alone is typically larger than a few meters. The LSP could decay fast enough 
when there are relatively large trilinear couplings, in which case a large trilinear contribution 
to a decay mode can dominate over a bilinear contribution making it impossible to measure 

One of the most promising ways of measuring through various L-violating channels is 
to look for same-charge dilepton events klj + 2W or lilj + 4j, from two LSP's produced 
in, e.g., e~^e~ collisions. As will be discussed in the following sections, the recent Super- 
Kamiokande data imply that one should observe comparable numbers of the same-sign (/i/x), 
(/ir), and (rr) events P] except for large tan/3 and small 



III. NEUTRINO OSCILLATIONS AND L VIOLATING PARAMETERS 

Let us estimate the sizes of the R-parity and L violating parameters for which the cur- 
rent neutrino data can be accommodated. In particular, our discussion is based on the 
assumption that R-parity violation gives rise to the neutrino masses and mixing explaining 
the atmospheric and solar neutrino oscillations. 

As is well known [^,TB-^, the mixing between neutrinos and neutralinos gives rise to 



the so-called tree mass. The neutrino-neutralino mass matrix Mn has two massless modes 
corresponding to the lightest neutrinos and the only massive neutrino which is clear from 
the following seesaw-reduced rank-one neutrino mass matrix; 

<^ = . (11) 

Let us first suppose that the tree mass gives rise to the largest neutrino mass m^g that is 
responsible for the atmospheric neutrino oscillation, namely, m^g = yAm^^^ ^ 0.05 eV P|. 
Then the overall size of ^ = Y,i it is determined from Eq. ( pT]) : 



e = 0.7x10-^-^ _ (12) 

^ cp \Mz) VO.OSeVy ^ ' 

Given the mass matrix (0), one can find the following expressions for the oscillation am- 
plitudes [0]: 



sin^ 29 



atm 


C2 C2 
,^2 S3 






atm 









sin2 2^^f = 4gg, (13) 

'-2 

'-)■ 

12' 
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The first quantity is measured to be sin 29^^ ~ 1 in the Super-Kamiokande which imphes 



^2 ~ ^3 0- The second one is measured in the CHOOZ experiment for i/g disappearance 
experiment |^ from which one finds ^f/^^ < 0.05 in the region Am^^j^ > 2 x 10~^ eV^. 
Therefore, as realized in Ref. 0, colhders may provide independent checks of the Super- 
Kamiokande and CHOOZ resuhs by measuring the ratios through the LSP decays 
discussed in the previous section. In other words, combining the results from the neutrino 
experiment and collider signatures, one could prove R-parity violation as a source of neutrino 
masses. For this to happen, of course, the LSP has to decay inside a detector. 

When the LSP is heavier than the W boson and thus the two body decays Xi ~^ I'^W'^ 
are allowed kinematically, the decay length is typically less than 10 cm [see the next section] . 
For the LSP mass below M^, the LSP can have only three body decay modes. The LSP 
decay through W, Z exchanges can make the decay length smaller than 1 m for the LSP 
mass close to Mw even ignoring the effect of the trilinear couplings. With larger trilinear 
couplings, the LSP decay through sfermion exchanges can be detectable in colliders. As will 
be seen in the next section, the contribution of the trilinear couplings with A, A' ~ 10~^ 
to the three body decay rates becomes comparable to the bilinear contribution in the case 
of common sfermion masses, ifif = 300 GeV. For the explicit calculation, we will take the 
universal trilinear couplings. 

The trilinear couplings of the order of 10~^ are too small to be influential in neutrino ex- 
periments. More interesting case is when the trilinear couplings are large enough to produce 
the neutrino masses and mixing explaining the solar neutrino deficit p2|. To estimate their 



sizes, let us recall that the trilinear couplings generate the one-loop neutrino mass through 
squark and slepton exchanges . Two neutrinos which remain massless at tree level become 
massive through the loop contribution. The loop mass takes the form, 

loop _ 3m|(A-K^ , ml{Ar + fitp) ^ 

^ij — o 2-2 ^i33^i33 o 2~2 ^*33^i33 , 

where Ai,^^ are trilinear soft parameters and are typical soft masses of sbottom and stau 
fields. Here we keep only Aj33,A-33 assuming that they give the largest contribution which 
is the case when the trilinear couplings follow the usual hierarchy of the Yukawa couplings, 
that is, the couplings involving the third generations are larger than the others. This would 
be a consequence of a flavor symmetry dictating the hierarchies among the Yukawa and 
trilinear couplings ||23| , |2^ . Still, having the arbitrariness of the L violating parameters 
.^i, Aj33 and A^33, the loop mass can be made larger than the tree mass. But in the context of 
minimal supergravity [Il6| , ^p5| or gauge mediated supersymmetry breaking |2^^8| where 



the misalignment 7^ is generated radiatively, the tree mass is typically larger than the 
loop mass. This is our motivation to assume that the largest mass from Eq. ([TT|) explains 
the atmospheric neutrino oscillation, and the second largest mass due to the loop correction 
accounts for the solar neutrino oscillation. Suppose now that A233 gives the second largest 
eigenvalue niu^, which will be the case with the model discussed below. Then, one needs to 



have m^, = -i/Am^, ^ 2 x 10 eV for the MSW solution of the solar neutrino deficit [122 



''1^2 ~ y "'"-sol 
Therefore, the size of A233 has to be 




5 / X i 



/i| / V2 X 10-3eV 
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taking = 300 GeV and ignoring the trilinear soft terms A. Let us note for the later use 
that the required value of A233 becomes smaller for larger tan/5. It is expected that A133 
is smaller than A233 by one or two orders of magnitude in the case of the favorable small 
mixing MSW solution As discussed in the previous section, one can use the decay 



modes Xi ~^ ^Ulj to get information on the flavor structure of Xijk- 



With such large trilinear couplings (p!5D , the three body decay rate can become even 
larger than the two body decay rate. To figure out when this can happen, let us note that 
the two body decay length for M^o ~ 150 GeV is about 0.1 cm [see FIG. 1]. Then using the 
Dawson formula for the three body decay rate of the photino-like LSP [^, we find that the 
same decay length can be obtained when a single large coupling takes the value, 

A.A'~lxlO-f^y7l5^W^y, (16) 

V300GeV/ \ M^o j \ I ) ' ^ ' 

where is the mass of squark and I = 1/T denotes the LSP decay length. We find the 
above value is in the range of the required value (|1^) for small tan j3. Therefore, the effect 
of trilinear couplings could surpass that of bilinear couplings so that it would be hard to 
separate out the information on alone. Furthermore, the loop masses from the trilinear 
couplings A^33 in Eq. (|T^) may affect the largest neutrino mass rriy^. For instance, taking the 
value of A^33 as in Eq. ([T3|), one finds that the loop mass is comparable to the tree mass ([TT]). 



Therefore, the expressions (|TT]), ([Tq) should be altered to include the loop mass as well. In 
the context of supergravity where nonaligned are generated from the trilinear couplings, 
the loop mass can be comparable to the tree mass to account for both the atmospheric and 



solar neutrino oscillations ^0[- this case, it would be nontrivial to probe the mixing angle 
through L violating decays of the LSP unless all the trilinear couplings are known. However, 
there is a class of models which are more predictive and do not possess the above difficulties. 
This is the so-called bilinear model [[r|, p!6| , p!7| , p6| ,p 



The bilinear model does not have trilinear R-parity violating couplings. A classic example 
would be the minimal SU(5) grand unified theory where generalized /i terms, fi^HuLa, are 
allowed but the trilinear couplings cannot be present to ensure proton longevity. As discussed 
in the previous section, in order to calculate the LSP decays, one has to go to the mass basis 
diagonalizing the neutralino/neutrino, chargino/charged lepton and slepton/Higgs mixing 
T6| . This amounts to generating effectively L violating trilinear couplings. Not bothering 



with details about bilinear and scalar mass parameters, it is expected that the trilinear 
couplings take the values, e.g., A^^^ ~ ^idjumj / {H^ . The dominant couplings are then 



a:33 ^ 2 X 10- 



8 1 / -TTV 



C$ \Mz 



As was discussed in the previous section, these couplings are negligible in our discussion for 
low tan (3 but become important for large tan/3. The other couplings are negligible for all 
tan P in view of our previous discussions. An interesting feature is that the loop mass can 
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account for the solar neutrino deficit through the MSW solution for large tan (3. For instance, 
taking tan/5 = 50, we get Aj33 ?s 2 x 10~^ which agrees with the value given in Eq. ([T5|) . We 
note that, in the large tan/3 region, the value of ^ in Eq. ( ]T2| ) becomes comparable to A233. 
In this case, the loop mass matrix has a large contribution from X'-^^ but it also takes the 
form oc ^j^j with the change of its overall size. Thus, the relations in Eq. (|13|) remain 
untouched. As will be shown explicitly in the next section, the induced trilinear couplings 
are large enough to make the LSP decay visible in colliders, but not large enough to change 
the two body decay length. Furthermore, the numbers C,i can be measured without any 
contamination from the trilinear couplings. 



IV. NUMERICAL ANALYSIS OF LSP DECAYS 

Before presenting our numerical results, let us recapitulate the crucial points for testing 
neutrino oscillations through the detection of the R-parity violating LSP decays in colliders. 
First of all, the LSP has to decay inside a detector. So, we require the LSP decay length / 
to be less than 1 m and show the corresponding region of neutralino mass parameters for 
given L violating parameters and sfermion masses. We also assume the unification relation 

Ml = ^4M2 , (18) 
and restrict ourselves to the region satisfying the recent LEP constraint on the mass of the 



lightest chargino: M~± > 90 GeV |31]. These restrictions rule out the region of M2 and fi 



for M^o < 50 GeV and the possibility of the chargino LSP. 

Since the mode hW or lijj for the two or three body decay case is important to determine 
the parameter [see Eqs. (i|,|ll])], we will present several figures showing the branching ratios 
of all possible decay modes on the M2 — /i plane, as well as with fixed M2 and varying /i. 
In the definition of e.g. the kjj branching ratio, we include both if. The quantities 
are measured by the bilinear contribution to the decay processes we need to separate out 
this from the trilinear contribution. For this purpose, we will present the branching ratios 
coming from bilinear terms and trilinear terms separately. To get an idea of when the 
trilinear contribution becomes visible, we will show the decay length and branching ratios 
taking the universal trilinear couplings of 10~®. As discussed in the previous section, when 
the trilinear couplings are of the order 10~^ which may be required to produce the solar 
neutrino mass, the trilinear contribution to the two body decay rate becomes comparable 
to the bilinear contribution. In this case, it would be meaningless to separate out the 
two contributions making it complicated to probe the neutrino oscillations through collider 
signals. However, it is possible to have a good separation between the bilinear and trilinear 
contributions while generating both the atmospheric and solar neutrino masses. This is the 
case in the bilinear model with large tan (3, where the order of the magnitudes of the induced 
trilinear couplings are given in Eq. (|1^. Therefore, taking the bilinear model as a specific 
example for the atmospheric and solar neutrino oscillations, we will show the decay length 
and branching ratios. 

In our computations, we take the following input values as a reference: 
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6«6 = e3 = ^, (19) 

A J = TJij = 300 GeV for all sfermions , 

while ^ fulfills the relation (|12|). We vary M2 and /i in the range of 50 ~ 500 GeV and 
—300 ~ 300 GeV. Here, we take ^2 = ^3 reflecting the large mixing angle of the atmospheric 
neutrino oscillation as in Eq. (13). 

Let us now present our results in FIGs. 1-8. In FIG. 1, we show the decay lengths 
(a) including only the LSP decays to VF, Z bosons or through Z exchanges, (b) taking 
the universal trilinear couplings A, A' = 10~^, and (c) considering the bilinear model, for 
tan/? = 3. The thick lines in the M2 — fi planes correspond to the points for M^o = M^y. 
The shaded region is excluded by the constraint on the lightest chargino mass. One finds 
that there is a very small region with the two body decay length larger than 1 m (denoted 
by filled stars just above the thick line) for negative /x. It is clear from FIG. la and Ic that 
the trilinear couplings in the bilinear model ( |T7|) are negligible for the small tan (3. There is 
also a little change in FIG. lb. 

The situation changes for tan/3 = 50 as in FIG. 2. The decay lengths become shorter 
than 10 cm in the two body decay region. The M2 — n plane is completely covered for the 
bilinear model as shown in FIG. 2c. This is mainly due to large A^33 as in Eq. (|T7|) . Also 
more region is covered in FIG. 2b with A, A' = 10"^. 

In FIG. 3, we present (a) the IW branching ratio adding all charged leptons, and (b) 
each branching ratio fixing M2 = 300 GeV, in the case of tan/3 = 3. We do not show the 
region with M^o < Mz for which the uZ channel is forbidden. As can be seen in FIG. 3a, 
there is a limited region (denoted by dots) where the IW mode is strongly suppressed, which 
corresponds to the point with fi = —140 GeV in FIG. 3b. This implies that a cancellation 
among various terms in 9^1^ in Eq. (|]) takes place. As we put ,^2 = (,3, the fiW (solid 
line) and tW (dashed line) branching ratios exactly coincide for large and thus for large 
M^o (M2 = 300 GeV). But for smaller \fi\ the two ratios deviate from each other which 
contradicts with the result of Ref. 0. This is a consequence of the charged lepton mass 
dependence of 9^^ in Eqs. (P) and (^ which becomes very important for the Higgsino-like 
LSP. Let us mention as a reference that the valley at ^ = —140 GeV corresponds to the 
LSP mass M^o = 123 GeV. Our result is larger than that of Ref. P in the uZ decay rate by 
factor of 2, coming from the fact that the AL = —1 and AL = 1 modes are to be summed. 
The eW branching ratio is negligibly small as we put ^1 << ^2,3- If Ci is comparable to 
^2,3, the curve for the eW branching ratio is in parallel with that of the nW branching ratio 
with the constant deviation by factor of ^1/^2 dependence on me or is negligible. 

However, the universal property does not remain for the eW and tW modes especially for 
large tan /3. 

The same set of the branching ratios as FIG. 3 are shown in FIG. 4 with tan/3 = 50. 
The IW branching fraction becomes larger than that with smaller tan /3 and it can be bigger 
than 40 %. From FIG. 4b, one finds that as tan/? increases the deviation between fiW and 
tW branching ratios greatly magnifies as discussed in section II. 

FIG. 5 shows three body branching ratios for Xi ~^ vvv., I'll, vjj and lijj without 
trilinear couplings. M2 is fixed to be 150 GeV. For tan /3 = 3 in FIG. 5a, the ^jj (solid line) 
and rjj (dashed line) branching ratios appear to be the same. On the contrary, there is a 
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big deviation between two branching ratios for tan (3 = 50 and small as shown in FIG 5b. 
This is also due to the large dependence of 9^^ on the lepton masses for the Higgsino-like 
LSP. 

In FIG. 6, three body branching ratios with universal trilinear couplings of 10~^ with 
M2 — 150 GeV and tan P — 3 are shown by separating (a) the bilinear and (b) the trilinear 
contributions. As there is almost no interference between the two contributions, the full 
branching ratios can be obtained just by adding these two contributions as in FIG 6c. The 
trilinear contribution to the Ijj modes is shown to be negligible for positive in FIG. 6b. As 
can be seen from FIG. 6d, the bilinear contribution (solid line) dominates over the trilinear 
contribution (dashed line) except for ~ —100 GeV. 

For tan /5 = 50 in FIG. 7, trilinear contribution becomes larger but still the bilinear 
contribution to lijj mode dominates. One again sees the big deviation between /ijj and rjj 
modes for smaller 

Finally, the three body branching ratios of the bihnear model with tan P — 50 are shown 
in FIG. 8. Here the trilinear contribution dominates but it comes mostly from the ujj mode. 
This is because the coupling A^gg do not contribute to the mode lijj since it involves heavy 
top production. Therefore, one can directly measure the lijj branching ratios from the 
bilinear contribution which have branching ratios of 1 — 6 % as shown in FIG. 8c. Another 
important property is that the vll branching ratio from the trilinear contribution is about 
1% and thus one can also probe the flavor structure of the trilinear couplings A. As a 
consequence, collider experiments can provide an independent check of whether the bilinear 
model explains the solar and atmospheric neutrino oscillations. 

V. CONCLUSIONS 

Recent observation of neutrino oscillation in Super-Kamiokande would require physics 
beyond the standard model. Supersymmetric extension of the standard model allowing L 
(and thus R-parity) violation is one of the well-motivated theoretical frameworks for nonzero 
neutrino masses and mixing. That is, the current atmospheric and solar neutrino data may 
be explained by R-parity violation in the MSSM. A distinctive feature of this picture is that 
it could be tested in colliders through L violating signatures of the LSP decay. In generating 
neutrino masses from R-parity violation, one has to take into account both trilinear and 
bilinear terms in the superpotential and in the soft breaking scalar potential. Typically the 
bilinear contribution to neutrino masses through nonaligned sneutrino VEVs dominates over 
the trilinear contribution. Based on this aspect, we have analyzed collider signals for the L 
violating decays of the LSP (taken as the lightest neutralino) which has two body and/or 
three body decay channels depending on the LSP mass. 

The results of our analysis are shown in FIGs. 1-8 presenting the decay lengths and 
branching ratios as functions of the supersymmetric parameters M2 and n under the as- 
sumption of gaugino mass unification. When the trilinear L violating couplings are too 
small, there is some region of the LSP mass where the LSP decay cannot be seen in a de- 
tector. However, if the trilinear couplings are large enough to give rise to the solar neutrino 
mass while the atmospheric neutrino mass comes from the bilinear parameters, the LSP is 
shown to decay inside a detector for the whole R-parity conserving MSSM parameter space 
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restricted by the recent lower bound on the chargino mass m^- > 90 GeV and the unifica- 
tion relation of gaugino masses. As a specific example of realizing both the atmospheric and 
solar neutrino oscillations, we have taken the bilinear model with large tan (3 which has no 
arbitrariness in the trilinear couplings. In such a model, the neutrino oscillation parameters 
measured in the atmospheric neutrino experiments can be directly probed in colliders by 
measuring the branching ratios of the LSP decays into kW or lijj. This opens a possibility 
to test the assumption of the R-parity violating MSSM as the origin of neutrino masses and 
mixing. 

In a more general case with arbitrary trilinear couplings, it would be impossible to 
directly test the neutrino oscillations in colliders. But if future coUiders measure the liW or 
lijj branching ratios it would provide a hint for or against a specific mechanism predicting 
the structures of the trilinear couplings and/or soft supersymmetry breaking parameters. 



VI. APPENDIX 



In order to express the amplitude of the decay x° 
introduce the notations: 



fifjfk in a compact form, we 



9f 



SwCW 



(20) 



for = —\ fermions , 



N*n^ ^^^^ for = i fermions , 



Pv{x) 

Qfix) = 



V2f 



2 X* 

w 



V2i 



oR 

w 



2^2 r2 



for V = Z,W, X = s,t,u , 



^ gJix-m^yrhjA^ _ 
\x — ml ]\x — rri; ) 

^ Jl J2 

[x — rrir )(x — ml ) 
_ {fn^f^ <-> m|^ ) ior X — u , 



Here m| r are the eigenvalues of the sfermion mass-squared matrix, 

where /S.j = mf{Af + iitp) or mf{Af + fi/t^) for = =i=i sfermions, respectively. 

For the amplitude of the decay x° — > j^ie'j'e^ coming from the bilinear contribution and 
the trilinear coupling Xijk, we find the following 5 nonvanishing matrix elements: 
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Qvee 
Qvee 



S^krfPw{t) 



LL 



- — A 



ijk 



PM-PeAu) 



oeud 
'-'LL 

oeud 



Pds) - IpuM - IPdii^) 



eud 
RL 



eud 
LL 



+^Kjk Pu,{t) - Paiiu) 



g ijk 

We find that the decay x° — >■ i^idjdk has 4 nonzero matrix elements: 



Qvdd \/ 

•^LL — \jk 



1 



1 



VEr = -S,knf,.Pz{s) - -KjkQ^t) , 



(21) 



Similarly, there are 5 nonvanishing matrix elements for the decay x° e^Ujdk'- 



(22) 



(23) 



For the decay 



nr^-^Kjk [Pd,it) - Pdti^)\ ■ 

UiUjUk, we find that there are only two nonvanishing matrix elements: 

^lT = Sjkul^Pzis) , (24) 
^lT = -<^.-fc^f.cPz(«) . 

Finally, the decay x'^ ~^ i^iT^ji^k has only one nonvanishing term: 



V^L = ^jkul^Pzis) + 5,knlMi) . 



(25) 
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FIGURE CAPTIONS 



FIG. 1. Two body and three body decay lengths of the LSP on the M2-H plane with 
(a) the vanishing trilinear couplings, (b) the universal trilinear couplings of A, A' = 10~^, 
and (c) the trilinear couplings given by the bilinear model when tan (3 — 3. The shaded 

area is excluded by the restriction m^± > 90 GcV and the two body decay is allowed only 
in the above the thick solid line corresponding to m^o = Myv- Different markers are used 
to denote the region with the decay length of / < 10~^ m (dot), 10^'^ < / < 10^^ m (open 
circle), 10~^ < / < 10~^ m (asterisk), 10~^ <l < \ m (filled circle), I > 1 m (filled star) and 
I <lm (filled box). 

FIG. 2. The same as FIG. 1, but with tan/? = 50. 

FIG. 3. (a) Two body branching ratio of the LSP decay into IW on the M2-11 plane 
with tan /3 = 3. Different markers are used to denote the region with the branching ratio of 
~ 0.2 (dot), 0.2 ~ 0.4 (open circle), 0.4 ~ 0.6 (asterisk), 0.6 ~ 0.8 (filled star), 0.8 ~ 0.9 
(filled triangle) and 0.9 ~ 1 (filled box), (b) The same branching ratios when M2 = 300 
GcV. In the shaded area, M^o < Alz- The four different lines correspond to the decay mode 
of Xi ~^ (solid line), Xi tW (dashed line), x5 ~^ (dash-dotted line) and ~^ 
(dotted line), respectively. 

FIG. 4. The same as FIG. 3, but with tan (5 = 50. 

FIG. 5. (a) Three body branching ratios with the vanishing trilinear couplings when 
M2 = 150 GeV and tan (3 — ?>. The shaded area is excluded due to the constraint M~± > 90 
GeV. The five hues correspond to the decay mode of Xi (sohd line), Xi — 'Tjj (dashed 

line), Xi i^jj (dash-dotted line), x? ^ ^^^^^^ (thick dash-dotted line), and x? — (dotted 
hue), (b) The same as (a), but with tan/? = 50. 

FIG. 6. Three body branching ratio with universal trilinear couplings of 10~^ when 
tan P — 3 and M2 — 150 GeV. Bilinear and trihnear contributions are shown in (a) and (b), 
respectively, and the sum of these two contributions is shown in (c). The line convention for 

each decay mode in (a), (b) and (c) is the same as in FIG. 5. In (d), we show the branching 
ratio summed over the five decay modes. The solid line denotes the bilinear contribution 
and the dashed line the trilinear contribution. 

FIG. 7. The same as FIG. 6, but with tan P = 50. 

FIG. 8. The same as FIG. 7, but with the trilinear couplings given by the bilinear model. 
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